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ABSTRACT 

Using a uniform analysis procedure, we measure spatially resolved weak gravita- 
tional lensing and hydrostatic X-ray masses for a sample of 18 clusters of galaxies. 
We find a radial trend in the X-ray to lensing mass ratio: at r25oo we obtain a ratio 
Mx/Ml = 1.03 ± 0.07 which decreases to Mx/Ml = 0.78 ± 0.09 at rgoo. This differ- 
ence is significant at 3(t once we account for correlations between the measurements. 
We show that correcting the lensing mass for excess correlated structure outside the 
virial radius slightly reduces, but does not eliminate this trend. An X-ray mass un- 
derestimate, perhaps due to nonthermal pressure support, can explain the residual 
trend. The trend is not correlated with the presence or absence of a cool core. We also 
examine the cluster gas fraction /gas and find no correlation with M^, an important 
result for techniques that aim to determine cosmological parameters using /gas- 
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1 INTRODUCTION 

Rich clusters of galaxies host the most massive collapsed 
dark matter halos. The X-ray emitting intracluster medium 
(ICM) bound to this halo can be a useful tracer of the dark 
matter content of the cluster. Under the assumption of hy- 
drostatic equilibrium the gradients of the gas pressure P 
and total gravitational potential $ are related by the simple 
differential equation 



VP = PgV$ 



(1) 



where pg is the gas density. The hydrostatic method is in 
principle powerful, because spatially resolved measurements 
of the gas pressure can constrain the shape of the halo and 
hence yield useful limits on fundamental dark matter prop- 
erties and cosmology. 

In reality, we do not know the reliability of hydrostatic 
mass estimates. Not only are merging clusters of galaxies in a 
nonhydrostatic state, but the plasma in many apparently re- 
laxed systems may be affected by additional nonequilibrium 
processes, which serve to boost P and hence cause an under- 
estimate of the cluster mass from X-ray observations of the 
thermal bremsstra hlung em i ssion. In early hydrodynamic N- 
body simulations, lEvrar 

3 (|l990l ) found the first signs of 
this underestimate, attributing it to incomplete thermal- 
ization of the ICM. More recent N-body work has shown 



that energy input from active galactic nuclei, pressure sup- 
port from turbulence and residual bulk motions, and varia- 
tions in the merging history may contribute substantially to 
this s ystematic bias (iDolag et aL 20051 : Faltenbacher et al 
'20051; iRasia et al.1 l2006l : iNagai et al.l l2007l : lAmegUo eTal 
200'g). Importantly, constraints on dark energy from the 
cluster mass function are highly sensitive to intrinsic scat- 
ter and systemat ic error in the mass-observable relation 
(I Lima fc "llull2005h : the scatter also significantly influences 
the normalization of the primordial fluctuation spectrum erg 
(Balogh et al.. 2006). 

Comparison with independent methods can provide a 
powerful means of checking the reliability of X-ray mass 
measurements. Here we focus on comparison of hydro- 
static X-ray masses with those derived from weak grav- 
itational lensing. Unlike the hydrostatic method, gravita- 
tional lensing does not require assumptions regarding the 
dynamical state of the cluster. Early work comparing lens- 
ing and X-ray masses showed that while strong gravita- 
tional lensing masses in merging clusters sometimes ex- 
ceeded X-ray masses by factors of « 1.5 — 2, weak lens- 
ing masses measured at larger radii are con sistent with 
the X-ray data to within » 20% l|Mira lda-Escu de fc Babul 



|1995| : ISauires et all 1 19961 : lAUen et al.1 il996. : .Squires et al 
199'^; AllenI 19981 ). This picture is supported by subse- 
quent ground- and space-based observations of relaxed 
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and merging clusters, albeit with lai r ly large uncertainties 
ll Allen et al.lT2 002: S mith et al.ll200ll: iHoekstra et al.ll2002l: 
Ettori fc Lombardil 120031 : ICvpriano et al.ll2004l : iHicks et all 



20061 ) 



Only recently has data with the accuracy required to 
carry out a systematic comparison of larger samples of hy- 
drostatic and lensing masses become available. Recent stud- 
ies, how ever, are a pparently in conflict. For a sample of 30 
clusters. [Pedersen fc Dahle (.2007, ) report a 30% excess in the 
normalization of lensing mass-temperat ure (M-T) relativ e to 
the X-ray value withiijj rsoo. However. iHoekstral (|2007| ) us- 
ing 20 clusters which also form the basis of this paper, report 
no significant excess in the M-T normalization at r25oo- 

To study the difference in X-ray and lensing masses 
in greater detail, we present a sample of 18 clusters — the 
largest yet with accurate spatially resolved hydrostatic as 
well as weak lensing masses. In our study, we take partic- 
ular care to account for possible systematic effects in the 
lensing and X-ray measurement process. Our data support 
a picture in which the lensing excess is negligible at r25oo, 
but increases at at larger radii. In Section [2] we discuss our 
data reduction procedure; in Section |3] we discuss the trend 
in the Mx — Ml relation; in Section |4] we account for po- 
tential biases; in Section [5] we consider trends in the gas 
fraction; and in Section |6] we conclude. We assume Ho = 70 
km s"^ Mpc"^ Q,M = 0.3, and = 0.7. 



2 DATA 

2.1 Lensing data 

The c lusters in our sample were drawn from IHoe kstral 
(|2007l ). which contains a weak lensing analysis of CFH12k 
data from the Canad a-France-Hawaii Telescope. We refer 
interested readers to IHoekstral (|2007t l for details of the 
data reduction and weak lensing analysis procedure. The 
sh ear measu r ement s discussed in this paper are identical 
to IHoekstral (|2007l ). but here we present updated values 
for t he masse s. The changes are due purely to the fact 
that IHoekstral (|2007l ) used the Hubble Deep Fie ld (HDF) 
redshift distribution ijFernandez-Soto et al.|[l999l ). which is 
based on data tak en within a much smaller field than the 
lllbert et al.] l|2006l ') study, based on the CFHT Legacy Sur- 
vey data. The conversion of the lensing signal into a physical 
mass estimate depends directly on the source redshift dis- 
tribution. This dependence is quantified by the parameter 
/3iens = max[0, -DjsZ-Ds], where Dig and Ds are the angular 
diameter distances between the lens and the source, and the 
observer and the source. We note that /3iens is an impor- 
tant parameter because it is degenerate with the projected 
cluster mass. 

We calculate the new /?iens value s using the photo- 



metric redshift distributions derived by lllbert et al.l l|2006l ). 
The resulting val ues are listed in Ta ble [T] The mean red- 



shift derived by Ilbert et al.l (|2006| ) is higher (also see 
iBeniamin et al.ll2007l ) resulting in a ~ 10% reduction in the 



inferre d cluster ma s s com pared to those in |Hoekstr3|2003). 
As the lllbert et al.l (|2006l ) results are based on four indepen- 
dent pointings, we can also examine the field-to-field vari- 
ation in /3icns. We find this to be a negligible effect, with a 
dispersion of only a 1—2 percent. Alt hough systematic b iases 
may still be present, the results from lllbert et al.] (|2006l ) pro- 
vide a significant improvement over previous studies. Conse- 
quently, we expect that the current uncertainty in the source 
redshift distribution afi'ects our mass estimates only at the 
few percent level. 

We base our lensing masses on the aperture mass es- 
timat es (for details see the discussion in §3.5 in IHoekstral 
IIoot'). This approach has the advantage that it is practically 
model independent. Additionally, as the mass estimate relies 
only on shear measurements at large radi i, cont amination 
by cluster members is minimal. IHoekstral |20o3) removed 
galaxies that lie on the cluster red-sequence and boosted 
the signal based on excess number counts of galaxies. As 
an extreme scenario we omitted those corrections and found 
that the lensing masses change by only a few percent at 
most. Hence our masses are robust against contamination 
by cluster members at the percent level. 

The weak lensing signal, however, only provides a direct 
estimate of the projected mass. To calculate 3D masses (such 
as M2500, MiQOQ, and M500) from the model-independent 2D 
aperture masses we project and renor malize a density profil e 
of the form ptot (r) (X r"^(r2oo + cr)"^ (|Navarro et al.lll997l '). 
The relationship between the concentration c and the virial 
mass is fixed at c oc M.^n'' ^^/(1 -I- z) from numerical simu- 
lations ijBullock et al.|[200ll ). Hence, the deprojection itself, 
though well motivated based on numerical simulations, is 
model dependent. 



2.2 X-ray data 

Data for 18 of the 20 clusters in Hoekstri diool) are available 
in the Chandra X-ray Observatory (CXO) public archive. 
We analyze and fit these data using the Joint Analysis 
of C luster Observations (JACO) package (|Mahdavi et al.l 
l2007f ). The reduction proced ure follows the detailed descrip- 
tion in lMahdavi et alT(|2007l ). except as discussed below. We 
briefiy summarize the procedure here. We reprocess the raw 
CXO data with CALDB 3.3, including the charge transfer 
inefficiency correction for ACIS-I data. To remove the par- 
ticle background, we subtract the appropriate blank-sky ob- 
servation renormalized to match 9-12 keV count rates in the 
most source-free region of the data. We examine the residual 
spectra in 9-12 keV energy range for each cluster to ensure 
that the particle background is cleanly subtracted. 

We then extract spectra in circular annuli centered on 
the X-ray surface brightness peak, masking any detected 
non-cluster sources. The residual diffuse astrophysical back- 
ground is fit as an additive component as described in the 
analysis below. 

The spectra are fit with a projected 3D cluster model. 
We begin with a mass profile consisting of an NFW dark 
matter distribution and a gas density of the form 



^ For a cluster at redshift 2, the overdensity radius is the ra- 
dius within which the mean matter density is A times the critical 
density of the universe at redshift 2. M/^ is the mass within that 
radius. 



-3ft /2 



(2) 



where rx,i are the core radii and /3i are the slopes of the Np 



© 0000 RAS, MNRAS 000, 000-000 



Mx/Ml relation 3 



Table 1. Weak-Lensing Measurements 



Cluster 




Miens 


-^2500 


'"2500 


-^1000 


^1000 




^500 


Abell 68 


0.255 


0.504 


2.56 + 0.60 


0.52 


4.51 + 1.30 


0.85 


6.64 + 2.44 


1.22 


Abell 209 


0.206 


0.627 


1.89 + 0.67 


0.48 


4.40 + 1.61 


0.86 


7.14 + 1.93 


1.27 


Abell 267 


0.230 


0.522 


2.14 + 0.39 


0.49 


3.70 + 1.10 


0.80 


6.29 + 1.93 


1.21 


Abell 370 


0.375 


0.412 


3.86 + 0.66 


0.57 


7.79 + 1.86 


0.97 


13.27 + 3.40 


1.47 


Abell 383 


0.187 


0.641 


0.76 + 0.33 


0.35 


2.29 + 1.04 


0.69 


4.49 + 1.96 


1.09 


Abell 963 


0.206 


0.592 


1.39 + 0.41 


0.43 


2.50 + 0.90 


0.71 


4.16 + 1.54 


1.06 


Abell 1689 


0.183 


0.639 


4.74 + 0.77 


0.65 


9.11 + 1.31 


1.10 


14.29 + 2.40 


1.61 


Abell 1763 


0.223 


0.574 


2.93 + 0.59 


0.55 


5.51 + 1.43 


0.92 


10.47 + 2.84 


1.43 


Abell 2218 


0.176 


0.644 


2.66 + 0.64 


0.54 


4.59 + 1.10 


0.88 


6.10 + 1.64 


1.22 


Abell 2219 


0.226 


0.562 


3.10 + 0.60 


0.56 


6.81 + 1.79 


0.98 


10.27 + 2.31 


1.42 


Abell 2390 


0.228 


0.586 


3.04 + 0.59 


0.55 


6.29 + 1.43 


0.96 


8.79 + 1.99 


1.35 


CL 0024.0+1652 


0.390 


0.379 


3.16 + 0.63 


0.53 


5.90 + 1.66 


0.88 


9.87 + 3.36 


1.32 


MS 0015.9+1609 


0.547 


0.267 


3.74 + 0.87 


0.53 


11.30 + 3.07 


1.03 


19.51 + 5.77 


1.56 


MS 0906.5+1110 


0.170 


0.674 


1.99 + 0.57 


0.49 


4.87+1.54 


0.90 


9.46 + 2.06 


1.41 


MS 1358.1+6245 


0.329 


0.447 


2.23 + 0.51 


0.48 


4.11 + 1.33 


0.80 


6.64 + 2.66 


1.18 


MS 1455.0+2232 


0.257 


0.564 


1.51 + 0.36 


0.43 


3.11 + 1.00 


0.75 


4.83 + 1.70 


1.09 


MS 1512.4+3647 


0.373 


0.434 


0.74 + 0.36 


0.33 


1.43 + 0.74 


0.55 


2.94 + 2.11 


0.89 


MS 1621.5+2640 


0.428 


0.368 


2.03 + 0.87 


0.45 


5.40 + 1.80 


0.84 


7.64 + 2.74 


1.19 



All masses are in units of IO^^Mq; all radii are in units of Mpc; z is the redshift of the clusters, 
and /3icns is a measure of the source redshift distribution. 



independent "/3-model gas distributions. We assume that the 
gas metallicity is of the form 

Z{r) = Zo{l+^/rl)-'"'^'\ (3) 

where Zq is the central metallicity, and rz and Pz describe 
the metallicity profile. 

The 3D temperature profile is calculated self- 
consistently using the equation of hydrostatic equilibrium: 

T{r) = - (^piooTioo + — dr j (4) 

where M is the total mass profile, /x is the mean molecular 
weight, nip is the proton mass, rioo is the radius at which 
the gas distribution is truncated, pioo is the gas density at 
the truncation radius, and Tioo is the temperature at that 
radius. We project the resulting 3D emissivity along the line 
of sight, absorb it by the gal actic hydrogen column as fixed 
bv lDickev fc LockmanI l|l99Cl ). convolve the result by the in- 
strumental response, and compare with the measured spec- 
tra using a statistic. A standard MEKAL plasma with 
variable abundance serves as the spectral model. 

The fitting procedure begins with a single /3 model fit 
(A^^ = 1); if this yields a poor fit, one or two /3 model 
components are added until the fit becomes acceptable. 
In two cases (Abell 1689 and Abell 2390) we could only 
obtain a good fit by also excising the central 70 kpc, as 
IVikhlinin et al.l l|2006l ) suggest. Other simultaneously fit pa- 
rameters are the gas metallicity profile (3 free parameters), 
the dark matter profile (2 free parameters), the normaliza- 
tions of the diffuse astrophysical backgrounds (3 free pa- 
rameters), and the temperature at the truncation radius 
Tioo- The astrophysical backgrounds are modeled as a power 
law with fixed slope 1.4 and free normalization (to model 
unresolved background AGN), and a unredshifted thermal 
plasma with free normalization and T < 0.5 keV (to model 
the local soft X-ray background). 

Accounting for the covariance of all measurable param- 



eters, proper treatment of the truncation of the gas distri- 
bution, simultaneous fitting of the background, a nd the ab- 
sence of subjectively weighted 2D temperatures (|Vikhlininl 
l2006f ) are unique features of the JACO code and are de- 
scribed in detail in iMahdavi et"al] (|2007l '). 

The measured masses are robust to variations of the gas 
truncation radius; varying this radius between r2oo and rso 
yields a 3% random variation in the measured masses. Tioo 
was often poorly constrained (only lower limits were possible 
with lower quality data). The JACO error analysis routines 
include this uncertainty in the final masses. 

The resulting total gravitating masses are listed in 
Table [2] In Figure [T] we compare these masses with the 
values reported for the same Chandra o bservations in 
iLaRoaue et"al] l|2006D and I Allen et~aLl (|2007l ). We find that 
our masses on the aver age in excellent agreem ent with the 
two samples, with the iLaRogue et alj 12009) values being 
somewhat lower, and the lAllen et al.l l|2007h values some- 
what higher, than ours. 



3 THE Mx /Ml RATIO 
3.1 A Radial Trend 

We examine the relationship between the X-ray and lensing 
masses at each density contrast using a simple constant of 
proportionality: 

Ma,x = a^M^^L (5) 

Both AIx and Ml are measured within the weak-lensing 
derived value of va- To estimate a a properly, we minimize 
a modified sta tistic appropriate ly weighted for errors in 
both coordinates (|Press et al.lll993) : 

2 {Ma,x — haMa.l)^ 

'^A,X + "'A'^A,L 
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Table 2. X-ray Measurements 



Cluster 


N0 






A-^iooo 






/2500 




Extrapolated 


Abell 68 


1 


98/98 


3.18 + 1.96 


4.98 + 4.10 


6.63 + 6.15 





09 ± 0.06 




1000,500 


Abell 209 


1 


184/179 


2.00 + 0.45 


3.97+1.77 


5.91 + 2.98 





12 + 0.03 




1000,500 


Abell 267 


1 


317/303 


2.27 + 0.31 


4.08 + 0.98 


6.26 + 2.05 





10 + 0.01 


0.030 + 0.030 


500 


Abell 370 


1 


915/804 


2.99 + 0.24 


5.81 + 1.47 


9.22 + 1.28 





12 + 0.01 




500 


Abell 383 


3 


402/366 


1.22 + 0.06 


3.04 + 0.65 


5.03 + 1.34 





10 + 0.01 


0.114 + 0.014 


500 


Abell 963 


3 


1129/1071 


1.92 + 0.21 


3.24 + 0.52 


4.20 + 1.64 





11 + 0.01 


0.079 ± 0.008 


1000,500 


Abell 1689 


3 


1065/917 


4.47 + 0.89 


7.40 + 2.32 


9.86 + 3.23 





12 + 0.02 


0.082 + 0.015 


500 


Abell 1763 


1 


440/419 


2.56 + 0.17 


4.78 + 0.50 


7.68 + 1.11 





12 + 0.01 




500 


Abell 2218 


2 


790/745 


2.43 + 0.49 


2.89 + 0.89 


5.99 + 2.33 





11 + 0.02 




1000,500 


Abell 2219 


1 


2272/2249 


4 62 + 47 


1 88 + 3 18 


1 7 32 + 6 50 





11 ± 0.01 




1 noo 500 


Abell 2390 


3 


3663/3341 


3.75 + 0.31 


7.53 + 1.09 


10.87 + 2.09 





14 + 0.01 


0.099 ± 0.009 


1000,500 


CL 0024.0+1652 


1 


165/127 


1.86 + 0.34 


4.87 + 0.74 


9.92 + 9.90 





09 + 0.02 


0.073 + 0.023 


1000,500 


MS 0015.9+1609 


1 


632/580 


2.67 + 0.19 


6.52 + 1.09 


10.28 + 2.42 





21 + 0.01 






MS 0906.5+1110 


2 


587/549 


1.59 + 0.19 


2.49 + 0.45 


3.23 + 0.73 





11 + 0.01 


0.030 ± 0.030 


500 


MS 1358.1+6245 


2 


752/675 


2.83 + 0.78 


6.31 + 2.81 


11.08 + 6.33 





09 + 0.02 


0.077 + 0.012 


500 


MS 1455.0+2232 


3 


1869/1685 


1.40 + 0.04 


2.25 + 0.10 


2.96 + 0.16 





15 + 0.01 


0.123 + 0.010 




MS 1512.4+3647 


1 


254/243 


0.69 + 0.13 


1.15 + 0.28 


1.74 + 0.44 





12 + 0.02 


0.095 + 0.028 




MS 1621.5+2640 


1 


237/255 


1.54 + 0.21 


3.61 + 0.54 


5.43 + 1.22 





13 + 0.02 







All masses at density contrast A are in units of IO^'^Mq and measured within the lensing radius rA as listed in Table [T] A^^ is 
the number of /3-models required to achieve a good fit; x^/*^ is the ratio of the statistic to the number of degrees of freedom 
u for the simultaneous fit to all spectra in all annuli; all masses are in units of 10^'* Mq. The cumulative gas fraction /2500 is the 
ratio of the gas mass to the total mass at r250o; ''cool is the radius (in units of Mpc) within which the isobaric cooling time is 
less than the age of the universe at the redshift of the cluster. Also shown are the density contrasts at which the masses had to 
be extrapolated beyond the available field-of-view. 



This formulation implies no intrinsic scatter in the data, 
an assumption we validate below. We derive errors in oa by 
locating the values at which X^^Xmin = li which correspond 
to the 68% confidence interval. 

The results appear in Figure [2l The simple one- 
parameter fit provides a good description of the relationship 
between Mx and Ml at all radii, and a clear trend is always 
present. The goodness-of-fit figure /v = 16/17, 11/17, and 
14/17 for masses measured within r25oo, r-iooo, and rsoo, re- 
spectively. Therefore, there is no evidence of intrinsic scatter 
in the Mx — Ml relation. All scatter can be explained by 
the statistical error. 

We find 02500 = 1-03 ± 0.07, aiooo = 0.90 ± 0.09, and 
1500 = 0.78 ± 0.09. The formal significance of the result is 
higher than it seems from the quoted errors, because the 
three slopes oa are correlated and therefore cannot be com- 
pared directly as a function of A. This is because the masses 
at each radius are correlated with the masses at other radii. 
Wherever the X-ray masses are extrapolated (as shown in 
Table [5}, they are highly correlated with M25oo,x; the lens- 
ing masses at all radii are highly correlated, because the 
aperture masses are derived by integrating over the tangen- 
tial shear from the radius of interest out to large radius. 

To properly evaluate the significance of the difference 
between a2500 and 0500, we employ a bootstrap procedure 
with, calculating 02500 — 0500 many times. The resulting dis- 
tribution appears in Figure |3] We find that 02500 — 0500 fol- 
lows a Gaussian distribution with a mean of 0.24 and a stan- 
dard deviation of 0.08. Thus, the difference in the Mx /Ml 
ratios at r25oo and rsoo is significant at the 2>a level. Us- 
ing a Monte Carlo simulation of the bootstrap processes, 
we confirm that the error obtained in this way is unbiased 
and the "error on the error" is close to the expected value, 
0.07/%/36 « 0.01, or 16%. 



3.2 Interpretation 

The variation of this constant of proportionality oa as a 
function of radius can yield useful constraints on cluster 
physics. This statistically significant gradient in Mx/Ml 
may be explained through an X-ray mass underestimate. 
Fo r example, the cosm ological N-body simulations discussed 
bv lNagai et al.l l|2007l ) predict a low Mx/Mtruc on account 
of nonthermal pressure support in the ICM: effects such as 
residual bulk motions cause hydrostatic mass estimates to 
be biased low. Earlier studies based on simulation s (e.g . 
l Evrardl Il990l: iFaltenbacher et al.1 l2005l : iRasia et al.) I2OO6I: 
ball man et al.1 120061 ) support this conclusion. iNagai et al.l 
(i2007l ) find that, on the average, clusters at 2 = ex- 
hibit little significant bias within r25oo (—0.12 + 0.16). How- 
ever, within rsoo, the fractional bias is more significant 
(—0.16 ± 0.10). Assuming that weak gravitational lensing 
mass has been corrected for projection effects (see below), 
the directio n and magnitude of the effect predicted by 
iNagai et all (|2007l ') at rsoo would seem to match our ob- 
served trend (see Figure |4]). 

An independent check on this result is a comparison 
with previously reported norma lization of the AIl ~ T re- 
lation. At r2500, iHoekstral l|2007h finds consistency between 
the normalizations of the Ml — T from CFH T data and the 
Mx -T rela tion from I Vikhhn in ct al.' ('2006'). On the other 
hand, at rsoo IPedersen fc Dahle (2007 ) finds the normaliza- 
tion to be high by « 30%. These results are consistent with 
the trend we detect here. 

If, as we suggest, incomplete thermalization of ICM is 
responsible for the X-ray mass underestimate, we should ex- 
pect clusters with the most substructure to exhi bit the great- 
est d eviations from Mx ~ Ml (e.g. Abell 781. ISehgal et all 
l2008h . Recent theoretical work has shown that it may be 
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Figure 1. Comparison of JACO measurements of M2500 witli 
other recent Cliandra measurement for a subset of 9 clusters. All 
masses are in units of W^Mp^ The Abell 383, 963, and 2390 
masses (filled points) are from Allen ct al. ('2007'); the remain- 
der (unfilled points) are from LaRoquc ct al. (2006). The JACO 
masses shown are estimated using X-ray data only, rather than 
calculated using the lensing values of r2500 as in Table [2] 



possi ble to correct for the X-ray mass underestimate in this 
way (jjeltema et al.ll2007l ). We look for such a trend in Fig- 
ure [5] by dividing the clusters into cool core and non-cool 
core subsamples, under the hypothesis that clusters with 
cool cores are less likely to have experienced recent mergers. 
We find that there is no statistically significant difference 
in the Mx — Ml relation for cool core and non-cool core 
clusters. This agrees with N-body simulations showing that 
the existence of a cool core is not a reliable pre dictor of the 
equil i brium state of a cluste r of galaxies (e.g. iPoole et al.l 
l2006l . l2007l : lBurns et allbOOTl ). More refined estimates of the 
level of substructure in each cluster are needed to evaluate 
the possibility of correcting for the X-ray underestimate. We 
will explore this avenue in a future paper. 
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Figure 2. Ratio of the hydrostatic X-ray mass to the gravita- 
tional lensing mass, both measured within the lensing overdensity 
radii r25oo, ''looOi a-^d rsoo. AH masses are in units of IO^^Mq. 
There is significant trend of decreasing Mx/M^ with increasing 
radius. Filled circles and unfilled squares show systems with and 
without cool cores, respectively (we define non-cool core clusters 
as those where the cooling time is nowhere smaller than the age 
of the universe at the redshift of the cluster). 
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Figure 3. The histogram shows the probability distribution of 
the difference between M^jMi^ measured at r2500 and Mx/Mi^ 
measured at rsoo- The difference is positive at 3(T significance. 
For comparison, we plot a normal distribution with ii = 0.24, 
o- = 0.08. 
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Figure 4. The Mx/Mi^ ratio as a function of overdensity A. 
The unfilled squares show the data before corre ction for excess 
struct ure along the line of sight as described by Ijohnston et al.l 
||2007|) : the filled circles show the data after correction for the 
effect. The trend is consistent with the X-ra y mass underestim ate 
predicted by N-body work (shaded region, iNagai et aIll2007^ . 



4 POTENTIAL BIASES 

As discussed above, we interpret the change in Mx /Ml with 
overdensity A as evidence of a deviation from hydrostatic 
equilibrium. However, an obvious concern is whether the 
result can be explained by biases in our measurements. Af- 
ter all, the geometries of real clusters are more complicated 
than assumed here (although the ensemble averaged system 
should be close to spherical). In this section we focus on a 
number of potential effects that could lead to a radial depen- 
dence. Note that other effects can change Mx/Ml, but do 
so independent of scale; for instance a change in the mean 
source redshift simply changes all Ma's by the same factor. 

4.1 Elongation Bias 

An explanation of our results could be that we are introduc- 
ing a systematic bias through spher ical modeling o f what 
we know ought to be triaxial systems. IPiffaretti et al.l (|2003| ) 
consider this question in detail. They find that for a broad 
distribution of elongations and inclinations, the differences 
between triaxial and spherical X-ray mass estimates are neg- 
ligible, of order 3%. This fraction is much smaller than the 
effect we observe. For clusters elongated nearly exactly along 
the line of sight, Piffarctti ct al. (2003.) do find that assum- 
ing spherical symmetry causes the projected X-ray mass to 
be underestimated by up to 30% at rsoo- It is conceivable 
that through X-ray selection, a number of clusters in our 
samples might be preferentially elongated close to the line 
of sight. However, even if present, s uch elongations c annot 
reproduce the trend we observe. As IPiffaretti et"al] (|2003l ) 
show, geometries that yield a 30% projected mass deficit 
at rsoo ought to show a similar deficit at r25oo, which we 
do not observe. Thus while the overall normalization of the 



trend in Figure |4] might be affected by preferential selection 
of elongated objects, the trend itself cannot be. 

4.2 Extrapolation Bias 

We also investigate whether the fact that most Msoo.x are 
extrapolated could affect the results. The four clusters for 
which no X-ray mass extrapolations are necessary are MS 
0016, MS 1455, MS 1512, and MS 1621. For these clusters 
alone, we find a25oo = 0.83±0.13, and asoo = 0.61±0.13. The 
difference in these ratios is similar to that of the full sample. 
Thus there is no indication that the substantial difference 
between 02500 and 0500 is caused by the extrapolation of the 
X-ray masses. For the extrapolated sample only, we obtain 
02500 = 1-07 ± 0.08 and 0500 = 0.84 ± 0.11. Thus the overaU 
ratios for the extrapolated and unextrapolated subsamples 
do not differ significantly (the subsamples are uncorrelated, 
and thus at fixed A we can fairly compare the nominal errors 
in OA for the extrapolated and unextrapolated subsamples) . 

4.3 Dependence on the Mass-Concentration 
Relation 

We check whether our use of the lBuUock et all l|200ll ) mass- 
concentration relation to deproject the weak lensing masses 
could cause the observed trend. We vary the normalization 
and redshift dependence of the mass-concentration relation. 
The corresponding mean concentration varies by ~ 30% 
(which we note is a rather extreme change), which leads 
to a systematic variation in the M25oo,l/M5qo,l ratio less 
than 5%. More recent N-body work shows that the slope 
of the mass-co ncentration rel a tion m ay be som ewhat shal- 
lower than the iBuUock et"aLl |200j) value (e.g. iNeto et al] 
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I2007I ). However, because our sample spans only an order of 
magnitude i n ma ss, the effect of changing the slope to the 
iNeto et al] (|2007l ') value (-0.1) is negligible. We therefore 
conclude that the use of an incorrect NFW concentration 
cannot reproduce a result of the strength we observe. 



4.4 Two-Halo Term Correction 

Finally, we consider the possibility that deprojections of 
weak lensing mass measurements in general are biased. 
When considering the accuracy of weak lensing masses at 
better than the 20% level, an important effect is the con- 
tribution of nearby large scale structure to the weak lens- 
ing mas^. The recent detection of a "two-halo term" — 
excess mass due to halos and uncoUapsed material outside 
the virial radius — in Sloan Digital Sky Survey (SDSS) clus- 
ters l|johnston et al.ll2007l ) provides direct observational ev- 
idence for considering such a contribution. 

The "two- halo" term could be an important effect in our 
measurements. We have carried out the deprojection of the 
weak lensing masses assu ming that they a re described by an 
NFW profile. However, in ljohnston et al.1 ^007) SDSS data, 
shear measurements out to 30 Mpc indicate that outside the 
virial radius, the matter density profile gr adually flattens. 
For th e range of masses we consider, the Ijohnston et al.l 
"two-halo" density profile is very approximately rep- 
resented by the form 



P2H oc 



l + (0.25r-/r2oo)' 
r{r2oo + crY 



(7) 



The above equation resembles an NFW profile within r2oo, 
but acquires a flatter shape outside 4r2oo. The profile is 
truncated at 30 Mpc, near wher e the cluster-mass corre la- 
tion function becomes neglig ible (|Havashi fc Whitell200'it l. 

If equation [7] describes the correct density profile, and 
yet we mistakenly use an NFW profile to deproject the lens- 
ing mass, then our 3D mass estimates would exceed the true 
value. It is straightforward to calculate exactly how much 
we would overestimate the masses; we can project both the 
NFW and the "two-halo" profile along the line of sight, and 
take the ratio of the estimated NFW and "two-halo" 3D 
masses. This ratio then represents a correction that we can 
apply to our data on a case-by-case basis. The results appear 
in Figure 131 We find that while the two-halo term correction 
does not remove the trend in aa, it does mitigate it, partic- 
ularly with regards to the large deficit seen at rsoo. 

The Mx/Ml trend remains because any extraneous 
mass along the line of sight affects mass estimates at all 
radii, and the difference between the excess at r25oo and 
rsoo is not sufficient to remove the radial trend in the un- 
corrected data. The corrected values are 02500 ~ 1.06 ±0.07, 
fliooo = 0.96 ± 0.09, and 0500 = 0.85 ± 0.10. If we repeat the 
bootstrap procedure in Section [3. II we find that the differ- 
ence between 02500 and asoo is reduced to 0.20 ± 0.08, still 
significant at 2.5(7, i.e. better than the 98% confidence level. 



^ As opposed to distant large scale structure which is a source 
of random noise and not a bias l lHoekstrall200l[) . 
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Figure 5. The X-ray gas fraction as a function of {top) weak 
lensing mass at r25oo and {bottom) redshift. There is no evidence 
of a trend with mass or redshift. Filled circles and unfilled squares 
show systems with and w ithout cool cores, r espectively. The solid 
horizontal line shows the lAUen et al value. 



5 GAS FRACTION TRENDS 

We also consider the X-ray gas fraction /gas as a function of 
lensing mass. The gas fraction is an important quantity for 
cosmology with clusters of galaxies. Under the assumption 
that /gas approaches a universal value at large enough radii, 
it is possible to d erive constraints on the dark energy equa- 
tion of state (see lAUen et al.ll2007l . and references therein). 
The constraints depend crucially on the assumption of uni- 
versality; if /gas were not a universal number independent 
of mass and re dshift, the dark energy constraints would not 
be correct (e.g. IVikhlinin et al.|[2006h . 

In Figure[5]we show the X-ray gas fraction as a function 
of the lensing mass at r25oo. There is remarkably little scat- 
ter in the relation. The clearest outlier in /gas is MS 0015.9, 
a previously known high /gas system with triaxial structure 
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elongated along the line of sight (jPiffaretti et al.|[2003l ). Oth- 
erwise, the gas fractions are closely clustered around 0.1. 
The mean value for all systems is {/2500) = 0.119 ± 0.006. 
This value i s in ex cellent agreement with the value found by 
lAllen et al.1 (|2007l ) for 42 clusters, /2500 = 0.110 ± 0.002. 

Importantly, /2500 is uncorrelated with both M25oo,i 
and z (Kendall's r is 0.079 and 0.217, respectively). We 
therefore cannot reject the hypothesis that /2500 is inde- 
pendent of lensing mass and redshift. 



6 CONCLUSION 

We provide a comparison of hydrostatic X-ray and weak 
lensing masses for a sample of 18 galaxy clusters. What 
makes our analysis unique is the uniform analysis of Chan- 
dra X-ray data and CFHT weak lensing data. This signifi- 
cantly improves the reliability of the comparison. At r25oo we 
find excellent agreement between lensing and X-ray masses: 
we obtain a ratio Mx/Ml = 1.03 ± 0.07. Interestingly, we 
observe a significant decrease in this ratio towards larger 
radii. At rsoo the ratio is Mx /Ml = 0.78 ± 0.09. Account- 
ing for correlations between the mass measurements, we find 
that the difference in the ratios is 0.24 ± 0.08, significant at 
3a. The trend of Mx/Ml with radius is consistent with 
previous measurements of the Ml — T relation, and is not 
caused by the assumption of spherical symmetry, extrapola- 
tion beyond the Chandra field of view, or uncertainty in the 
mass-concentration relation. 

The trend remains even after we correct for a systematic 
overestimate of the weak lensing mass due to correlated large 
scale structure. We show that the trend is consistent with 
simulations in which nonthermal pressure support causes a 
systematic underestimate of the cluster mass. Interestingly, 
the underestimate is not correlated with the presence or ab- 
sence of a cool core. These results are relevant for upcoming 
surveys that aim to measure the cosmological parameters 
from large cluster samples, especially for observations seek- 
ing measurements closer to rsoo. 

We also determine the X-ray gas fraction at r25oo; this 
quantity has been used in several studies to measure w, the 
parameter of the dark energy equation of state. We find that 
/2500 is not correlated with lensing mass or redshift. 
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